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The high-transition-temperature (high-Tc) superconductivity discovered recently in iron 
pnictides is analyzed within a perturbation theory. Specifically, the probable pairing sym- 
metry, the doping dependence of the transition temperature and the pairing mechanism are 
studied by solving the Eliashberg equation for multi-band (2- and 5-band) Hubbard models 
with realistic electronic structures. The effective pairing interaction is expanded perturba- 
tively in the on-site Coulomb integrals up to third order. Our perturbative weak-coupling 
approach shows that sufficiently large eigenvalues of the Eliashberg equation are obtained 
to explain the actual high transition temperatures by taking realistic on-site Coulomb inte- 
grals in the 5-band model. Thus, unconventional (non-phonon-mediated) superconductivity 
is highly likely to be realized. The superconducting order parameter does not change its 
sign on the Fermi surfaces, but it does change between the electron and hole Fermi surfaces. 
Consequently, the probable pairing symmetry is always "a nodeless extended s-wave sym- 
metry (more specifically, an s-t-wave symmetry)" over the whole parameter region that we 
investigated. It is suggested that the 2-band model is insufficient to explain the high values 
of Tc. 

KEYWORDS: iron-pnictide superconductors, superconducting mechanism, pairing symmetry, 
perturbation theory, multi-band Hubbard model 



1. Introduction 

The recent discovery of high-Tc superconductivity in iron pnictides and related compounds 
has strongly intensified the research activity in solid-state physics. The highest transition 
temperature attained in iron pnictides is about 56 K to date, which is currently the highest 
reported superconducting transition temperature among noncuprate superconductors. 

A prototypical iron-pnictide superconductor LaFeAsOi-xF^^ (x =0.05-0.12, Tc = 26 K) 
was discovered by Kamihara et al.^-* Before the discovery of the arsenic superconducting 
system LaFeAsOi-xF^;, it had been discovered that LaFePOi-^Fa;, where As in the arsenic 
system is replaced by P, is also a superconductor with Tc =3.2-6 K.^) The LaFeAsOi-^Fa; 
crystals are formed by stacking LaOF and FeAs layers alternately. Each iron atom is sur- 
rounded tetrahedrally by four As atoms, and each FeAs layer is constructed by connecting 
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the edges of the As tetrahedra. In the undopcd case x = 0, LaFeAsO appears to be a poor 
metal, and undergoes a structural transition from tetragonal to orthorhombic at 155 K and 
spin-density-wave (SDW) ordering below 137 K (with a wave vector q = (vr, vr, vr) and a mag- 
netic moment about 0.3-0.4/iB on each iron atom at K)."^) By substituting oxygen with 
fluorine, doped electrons cause metallic conduction with the disappearance of both the SDW 
and structural transitions for x > 0.05, and superconductivity occurs below about 26 K.^^ 
Moreover, subsequent studies clarified that replacement of La by other rare earth elements 
increases the transition temperature markedly: 41 K for Cc,^^ 52 K for Pr,^^ 52 K for Nd,*^-* and 
55 K for Sm7'^) Superconductivity has also been discovered in fluorine- free systems, includ- 
ing KFeAsOi_j/ (R=La, Ce, Pr, Nd, Sm, Gd. = 31-55 K),^-^^) Gdi_^Th^FeAsO (Tc =56 
K).i3) Hole-doped Lai -a^Sr^iFeAsO also shows superconductivity with Tc =25 K.^^^ Thus, a 
superconducting phase exists not only in the electron-doped side but also in the hole-doped 
side, and it appears that the superconducting transition is inhibited due to the structural 
distortion and SDW ordering only near the undoped case a; = 0. It is natural to expect that 
superconductivity in all of these compounds occurs basically by the same mechanism, in spite 
of the quantitative difference in Tc. 

In addition to J?FeAsOi_^F^ (R: rare-earth element) series, many other related supercon- 
ducting iron compounds without oxygen have been discovered, e.g., Bai_2;K2:Fe2As2 {Tc =38 
K),i5) Lii_^FeAs {Tc = 18 K), 16-17) ^nd Cai_^Na^Fe2As2 (Tc =20 K)}^^ In these related 
compounds, FcAs layers are separated not by layers of rare-earth elements and oxygen, but 
by layers of BaK, Li, or other various elements. Therefore, we naturally consider that the 
feature essential for the superconductivity is the existence of FeAs layers, while the remaining 
elements may play only a relatively minor role of a carrier reservoir or spacing between the 
FeAs layers. In fact, similar superconductivity has also been observed in a-FeSe (Tc =8 K),^^) 
in which the valence number of the iron atoms may be the same as that in iron pnictides but 
no spacing elements are inserted between the FeSe layers. This fact suggests that only the 
layers of square networks constructed from iron atoms may be necessary for superconductivity. 

Electronic structure calculations for some of the above-mentioned iron-based superconduc- 
tors have already been carried out in many studies. ^^^^^^ According to these studies, the elec- 
tronic states near the Fermi level are formed dominantly by the FeSd electronic states. Thus, 
iron pnictides may be considered to be typical strongly correlated electron systems, in which 
the strong Coulomb interaction among the iron 3d electrons can affect the electronic proper- 
ties significantly. The Fermi surface consists of hole pockets around the F point and electron 
pockets around the M point, qualitatively consistent with the results of angle-resolved pho- 
toemission (ARPES) experiments^^"^^) and more recent quantum oscillation experiments.^^) 

An important step toward understanding this newly discovered iron-pnictide supercon- 
ductivity will be to clarify whether it is conventional phonon-mediated superconductivity or 
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unconventional superconductivity such as that in copper oxides. There are already several 
reasons why unconventional pairing may be realized in iron pnictides: (i) Tc is very high, com- 
pared with conventional phonon-mediated BCS superconductors, (ii) Fe3d states dominate the 
majority of the density of states near the Fermi level, as mentioned above, and therefore 
the electron correlation originating from the Coulomb interaction on the relatively localized 
Fe3d bands is considered to be strong in these systems, and (iii) electron-phonon coupling is 
expected to be weak according to first-principles calculations.^^) Thus, it will be a very inter- 
esting issue to determine whether or not we can understand iron-pnictide superconductivity 
as a natural result of electron correlations, as in the case of unconventional superconductivity 
in many other strongly correlated systems. 

So many theoretical works on iron-pnictide superconductivity have already been presented 
that we cannot review all of them here. Many of them focus on the pairing mechanism and 
pairing symmetry. Among the proposed pairing symmetries, the most promising pairing sym- 
metry is a spin-singlet extended s-wave (s±-wavc) symmetry, in which the superconducting 
order parameter changes its sign between the Fermi surfaces. This pairing symmetry has 
been suggested in considerations based on the magnetic susceptibility,^^) random-phase ap- 
proximation (RPA),^^'^^) renormalization group (RG) calculations,^^' fluctuation-exchange 
approximation (FLEX)'^^) and other studies. The s±-wave pairing symmetry appears now 
to be consistent with the results of most experiments, including those based on ARPES,^''''^^) 
penetration depth, ^'^^ and neutron scattering. A recent theoretical study suggests that the 
sign reversal of the ,s±-wave pairing symmetry is consistent also with nonmagnetic impurity 
effects. Further investigations are still desirable to confirm that this pairing symmetry is 
indeed realized in iron-pnictide superconductors. Our present theoretical study will provide 
additional support for this pairing symmetry. 

In the present work we present a microscopic theory of iron-pnictide superconductivity. We 
take a weak-coupling approach on the basis of third-order perturbation theory with respect to 
the electron Coulomb interaction among the Fe3d electrons. This perturbative approach has 
been applied to many realistic electronic structures in order to discuss the pairing symmetry 
and pairing mechanism in real strongly correlated electron systems, and has succeeded in 
explaining various exotic pairing states, including the dy.2 _y2-wscve states in cuprate and or- 
ganic superconductors,^'^"^^) the triplet p-wave state in Sr2Ru04,^^'''^) the dr^2_y2-w&Ye states 
in several heavy-fermion super conductor s.^^~^^) Therefore, it is interesting to also apply this 
approach to iron-pnictide superconductors. Specifically, the probable pairing symmetry, the 
momentum and band dependences of the superconducting order parameter, the doping depen- 
dence of the transition temperature, and the pairing mechanism are investigated by solving 
the Eliashberg equation for 2- and 5-band Hubbard models with realistic electronic struc- 
tures. The effective pairing interaction is expanded perturbatively in the on-site Coulomb 
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integrals up to third order. Our weak-eoupling approach shows that sufficiently high values 
of Tc are obtained for realistic values of the on-site Coulomb interaction in the 5-band mod- 
els in order to explain the experimentally observed high values of Tc. Thus, it is concluded 
that unconventional (non-phonon-mediatcd) pairing is highly likely to be realized in iron- 
pnictide superconductivity. The superconducting order parameter does not change its sign 
on the Fermi surfaces without any nodes, but it does change between the electron and hole 
Fermi pockets. Consequently, the probable pairing symmetry is a spin-singlet extended s-wave 
or, more specifically, an s-^-wave symmetry. The calculated transition temperatures do not 
depend so sensitively on doping concentration, consistent with the results of experiments, but 
should drop suddenly around doping levels where some electron or hole pockets disappear. 
This suggests that nesting between the electron and hole Fermi pockets is responsible for this 
unconventional pairing. It is suggested that the 2-band model is insufficient to explain the 
high transition temperatures. The present article is a detailed version of our previous short 
article.^®) 

2. Formulation 

In the present work, we use multi-band tight-binding models to describe the electronic 
structure near the Fermi level. Taking account of the two-dimensionality of the actual crystal 
and electronic structure, we use two-dimensional tight-binding models. The unit cell in a single 
FeAs layer contains two iron atoms. Here we should note that the electronic structure near 
the Fermi level is dominated by the Fe3d orbitals, as many electronic structure calculations 
demonstrate, and will be reproduced effectively by taking only the Fe3d-like local orbitals (or 
more specifically, maximally localized Wannier orbitals, as presented in ref. 49). As Kuroki 
et al. pointed out,^^^ by extracting only Fe3d-like states, we can effectively construct a tight- 
binding model on a simple two-dimensional square lattice consisting of only iron sites. This 
simplified scheme helps us greatly by reducing the number of energy bands to be considered. 
Thus, the noninter acting part of the Hamiltonian is expressed as follows: 

i,j e,e' o- 

where i and j denote iron sites on the square lattice, i{£') denote FeSd-like localized orbitals, 
and Ciirr (4^^) is the annihilation (creation) operator of electrons of FeSd-like orbital i at site 
i with spin a. Here note again that the index specifying each of the two iron atoms in a unit 
cell is not necessary. Specific choices of tifjc are given in §3.1 and §4.1. 
We can diagonalize Hq by using a unitary transformation in the form 

^0 = EEE^«(^)4aCfcaa, (2) 
k a T 

where the index a specifies each of the diagonalized bands and k is momentum. Ea{k) is the 
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diagonalized-band energy, and Ckaa is related to cuo- by 

CUa = 5" 5" Uia{k)e'^-'''Ckaa, (3) 

where N is the number of iron lattice sites and Uia{k) represents the elements of the diago- 
nalization matrix. 

The electrons on the local orbitals will be affected by the on-site Coulomb interaction. 
The interacting part is given in the following form: 

J J' 

+ 2 Z Z Z 4a4v'Ci€<T'Ci^v + y Z Z Z cLcL'Ci€V'Ci€v, (4) 

where U ^ U', J, and J' are the Coulomb integrals for the intraorbital Coulomb interaction, 
the interorbital Coulomb interaction, the Hund's coupling, and the interorbital pair transfers, 
respectively. The total Hamiltonian is given by 

H = Ho + H'. (5) 

Propagation of the electrons of band a in the normal state is described by the normal 
Green's function: 

where A; = (fe, ia;„) and Un is the fermionic Matsubara frequency: a;„ = (2n + 1)'kT. Super- 
conducting transition points can be determined by the linearized Eliashberg equation: 

A„.,.,(fc) = -^ZZ Z yaa,a,,a'a,.,{k,k')\Gf{k'tl^a'a^a,{k'), (7) 



where ai are spin indices, T is temperature, Aa(j-^(j,^{k) is the anomalous self-energy on band 
a, and Vacricr2,a'a3(TA{kTk') is the effective pairing interaction. We have used the shorthand 
notation = J^feX^^^^- At superconducting transition points, the Eliashberg equation has 
a non-self-evident solution, i.e., /S.aaa'{k) ^ 0. Replacing the left-hand side of eq. (7) by 
\l^acn(T2{k) 1 then eq. (7) is regarded as an eigenvalue equation with the eigenvalue A and 
eigenfunction Aa^i (72(^)5 transition points are determined by the condition Amax = 1, 
where Amax is the maximum eigenvalue. Since Amax is usually a monotonically decreasing 
function of T, therefore the superconducting instability is already attained for Amax > 1 and 
normal state is still favored for Amax < 1- The momentum dependence and spin symmetry of 
the eigenfunction /^aa-^a2 (k) giving Amax characterize the most favorable pairing state. 

Evaluation of the effective pairing interaction V^o-i fT2,a' 0-3(74 (^j k') is a quantum many-body 
problem. The pairing interaction is expressed in the orbital representation by using the diag- 
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onalization matrix elements Uia{k) as follows: 

{in} 

where ("„ is the orbital-spin combined index, i.e., = (^tm ^n)i U^^{k) is the complex conjugate 
of U^aik), and the summation with {in} denotes that with respect to all of ii, £2, £3, and £4. 
Here we define the Green's function in the orbital representation by 

Gfiik) = E Ue,a{k)UUk)G^^Hk)Sa,a,. (9) 

a 

We note that H' can be expressed in the form 

-f^' = 9 XI X] -^ClC2,C3C44€lO-i4£20-2<^*^30-3Ci€40-41 (10) 



2 

i {Cn} 



with 



h 






= u 


h 






= U'-J 


h 




{Vci){lc) 


= U' 


h 






= J 


h 




e'a){e'a) 


= J' 



(11) 



and the other elements are zero, where £ £' and a a. The summation with {(n} denotes 
that with respect to all of the orbital and spin indices, in and cr„. It is useful to define the 
antisymmetrized bare Coulomb vertex by 

^Cikc3C4 = ^ClC2,C3C4 - ^ClC2,C4C3- (12) 

We expand V^if2,C3C4(^' perturbatively to third order in H', using G^^^^^{k) and T^'^^^ C3C4' 
Diagrammatic representations of the perturbation terms are displayed in Fig. 1. The Eliash- 
berg equation in the orbital representation is written using the anomalous Green's function 
^<iC2(^) as 

AciC2(^) = -^EE^C3C4(^')%C2,C4C3(^>^')- (13) 

k' C3C4 

Through the expansion of V"^iC2,C4C3(^' ^')' anomalous self-energy A^^^2(^) is also expanded 
as 

+Af-^\k) + Ag-f + Af-^'\k) + Ag-f (14) 
The analytic expressions of the contributions are 

4iW = -^§EE^C3C4(fc')rgk,,,3, (15) 

k' C3C4 
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(3-3) 



(3-3*) 





Fig. 1. Diagrammatic representation of the effective pairing interaction V^i {2,^3^4(^1 '^^ third- 
order perturbation theory. Each vertex and sohd line with an arrowhead represent the on-site 



Coulomb interaction T 



(0) 

ClC2,C3C4 



and the Green's function G^^^^^{k), respectively. The contributions 
enclosed by the broken rectangle, (1), (2), and (3-1), arc also included in RPA calculations, while 
the remaining contributions, (3-2), (3-2*), (3-3), and (3-3*) are the vertex corrections and are not 
included within the RPA. 



]^ E E E ^c»c.(*')x!,°!„,,m(* - *')r5' 

174,72C3^7lC2,C4T3 ' 



(16) 



k' C3C4{7n} 



N 



k' C3C4{7n}{en} 



Cl74,72C3 7l?4,673 ClC2,C4C3' 



(17) 



^^ / EEE E^C3C4(^')G(l(A.-fc' + A:,)xgUe3(^ + ^i) 

^ fc'fclC3C4{7n}{Cn} 
XLr^g^^^^/Clji ^i^4,-^2C3^ 7ie3,C4€2^ €^2,^473' 

- (^) ' E E E E ^C3C4 {k')G^L i-k + k' + h)xfi,,,s-k + h 

^ ^ k'ki C3C4 {7n} 
XLr^^^gi^ftlJi ^,^2,74«3^ ^472,aC3^ 73C2,C47l' 



(18) 



(19) 



1 /r 



-K^ EEE E^C3C4(fc')Ga,(-^+fc'+fci)4;U,3(-fc+^i) 

^ ^ fc'fcl C3C4 {7«} {€n} 
~»^7473l«;iJJ- Cl72,74C3^ 6e2,C471 73C2,a«4' ^"^^ > 
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^ ^ k'kl (3(4 {7n} {Cn} 
>^Ll^^^^{Kl)L ^j^^^^g^jJ- ^i^4,72C3^ 7lC2,C4T3' ^■^■^^ 



with 

x£k?3«.(^) = E 41(^)4S3(^ + k), (22) 

<k«3..(«) = E Gfi{k)GfXil - k). (23) 

The factor of ^ in eqs. (15), (20), and (21) is necessary to avoid double counting. The band 
indices are replaced with orbital indices by using the diagonalization matrix: 

^aa,a,{k) = J] t/;^„(fc)t/,;„(-fc)A^,^, (fc), (24) 

^aUk) =^Ue,a{k)Uhai-k)Faa^aM- (25) 
a 

The anomalous Green's function is related to the anomalous self-energy by 

Fa.i.2(fc) = |G'f (fc)l'A,.,.,(fc) (26) 

in the band representation. Note that this relation is valid only in the region where the 
superconducting order parameter is nonzero but small (i.e., in the vicinity of transition points), 
since all quantities in the present formulation are linearized with respect to the anomalous 
functions. It is clear that eq. (7) can be consistently obtained from eqs. (13), (25), (26), (24), 
and (8). The anomalous functions Aao-icr2(^) ^iiid Facna-zik) can be determined by eqs. (7) and 
(26), which is carried out numerically. For numerical calculations in the present study, we take 
32 X 32 fc-points in the first Brillouin zone and 512 Matsubara frequencies. The momentum 
and frequency summations are performed using the fast Fourier transformation. 

We now turn our attention to the remaining degrees of freedom in the phase factor of the 
diagonalization matrix elements Uia{k). We can easily verify that if Uiaik) diagonalizes Hq, 
then Uia{k)e^'^'^^''^ also diagonalizes Hq, where 4>a{k) is an arbitrary real function of k. It is 
easy to see that G^^^i^^{k) and V(^-^c,2,C,3C,i{k-,k') do not depend on the choice of phase function 
(j)a{k). However, the values of V^o-icr2,a'o-3cr4(fc, ^0 and the anomalous functions depend on 
(j)a{k) in general, since they are subject to extrinsic phase modulations. In order to avoid this 
unfavorable situation, it is natural to assume the following condition: 

U,a{k) = UU-k). (27) 

We use this condition throughout the present work. 

In the present formulation, we have taken only the diagonal components of anomalous 
functions with respect to band indices and have neglected the off-diagonal components. Here 
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note that the diagonal eomponents describe intraband pairing, while the ofF-diagonal ones 
describe interband pairing. Superconducting instability is generally connected with the loga- 
rithmic increase in the pair-correlation function at low temperatures. While the diagonal pair 
correlations exhibit this increase down to zero temperature, the off-diagonal ones do not in 
general. This is because the increase of the pair correlations between different bands stops 
below a temperature comparable to the energy difference between the bands. Thus, we may 
neglect the off-diagonal components for the reliable prediction of pairing symmetry at low 
temperatures. 

The same formulation presented here was previously applied to the spin-triplet supercon- 
ductivity in Sr2Ru04 by the present author and Yamada.^^^ In the present section we have 
given improved and more useful analytic expressions for the perturbation terms. 

3. Analysis of 5-band Hubbard model 

3.1 5-band tight-binding model by Kuroki et al. 

In this section, we present an analysis based on a 5-band model. Following Kuroki et al,^^) 
we take the five local Fe3d-like orbitals: £ = 3Z^ - B!^, XZ, YZ, -Y'^, and XY. The 
tight-binding parameters ti£j£i are provided in ref. 23. Here we should bear in mind that the 
original unit cell contains two iron atoms. The principal axes X and Y are taken along the 
basal plane parallel to the FeAs layers, and the Z-axis is taken to be perpendicular to the 
basal plane, with respect to this original unit cell. The above orbital states are represented 
using these original coordinates. Neglecting the As sites in the effective models, the remaining 
iron atoms form a two-dimensional square lattice in which each unit cell contains only one 
iron atom. The principal axes x and y with respect to this new unit cell are parallel to the 
nearest-neighbor iron-iron bonds and are obtained by rotating the original axes X and Y by 
45° in the basal plane. In momentum space, the original Brillouin zone in which positions 
are parametrized hy k = {kx, ky) is extended to a wide Brillouin zone in which positions are 
parametrized by A; = {kx, ky). 

The electronic structure of the 5-band model is shown in Fig. 2. The total electron number 
n is related to the doping level x hy x = n — 6. The Fermi surface for a; = 0.1 forms five Fermi 
pockets in the extended Brillouin zone: two hole pockets {a and /?) around (0,0), one hole 
pocket {a') around (7r,7r), and electron pockets (7 and 7') around (7r,0) and (0,7r). In the 
original folded representation, both (0, 0) and (tt, tt) in the extended Brillouin zone are folded 
onto the T point, and (tt, 0) and (0,7r) are folded onto the M point. Thus, three concentric 
hole-like surfaces will be located around the F point, and two concentric electron-like surfaces 
will be constructed around the M point, in the original folded Brillouin zone. The XZ, YZ, 
and X^ — states occupy most of the density of states near the Fermi level and are naturally 
expected to dominate the electronic properties including the superconductivity. 
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Fig. 2. (Color online) Electronic structure of the 5-band model by Kuroki et al. in ref. 23. (a) Fermi 
surface and (b) energy band dispersion, both of which are depicted in the extended Brillouin zone, 
(c) total density of states and partial density of states for the orbitals, 3Z^ — R^, XZ, YZ, X^ — Y"^ 
and XY . The Fermi level is set to zero on the horizontal axis. 



3.2 Favorable pairing symmetry 

The most favorable pairing symmetry is obtained from the eigenfunction giving the max- 
imum eigenvalues of the Eliashberg equation. In Fig. 3, the eigenvalues for various pairing 
symmetries at x = 0.10 are shown as a function of temperature. The most favorable pairing 
symmetry is the spin-singlet s-wave (or more appropriately, as we shall see below, the ex- 
tended s-wave or s±-wave), and the maximum eigenvalue is sufficiently large to explain the 
actual high transition temperatures for realistic values of the Coulomb interaction: U = 1.2 
eV, U' = 0.9 eV, and J = J' = 0.15 eV. In Fig. 3, the transition temperature is evaluated to 
be about 100 K from the condition Amax(^c) = 1- 

We number the energy bands Ea{k) (a = 1, • • • , 5) in ascending order with respect to 
energy, i.e., a < a' for Ea{k) < Ea'{k) at any k. The three bands, a = 2,3, and 4, cross 
the Fermi level. The momentum dependences of the anomalous Green's functions Fa{k) at 
Un = vrT on these three bands are displayed in Figs. 4(a)-(c). Here we have suppressed the 
spin part of the spin-singlet anomalous Green's functions by using the Pauli spin matrix 

<7j/: Faaiaiik) = Fa{k)[ 

i'7j/](TiCT2 • As we see in Fig. 4, the zero contours of the anomalous 
functions do not have any cross sections with the Fermi surfaces. Thus, we conclude that the 
superconducting gap will fully open on the Fermi surfaces, noting that the nodal positions 
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Fig. 3. (Color online) Eigenvalues of the Eliashberg equation as a function of temperature for various 
pairing symmetries. The doping level is x = 0.10. 



of the superconducting gap generally coincide with those of the anomalous self-energy. A 
remarkable feature is the sign change of the superconducting order parameter between the 
hole pockets (a, a' and /3) and the electron pockets (7 and 7'). Thus, the pairing symmetry 
is not a simple s-wave, but an extended s-wave symmetry, or more specifically, an s±-wave 
symmetry. 

The superconducting order parameter on the inner hole pockets a and a' around (0, 0) 
takes nearly the same absolute value as that on the electron pockets 7 and 7' around (vr, 0) 
and (0,7r), although the sign of the order parameter is changed between them. The absolute 
magnitude of the order parameter on the outer hole pocket /3 is about half of that on the 
inner hole pockets a and a'. These features are in agreement with the results of ARPES 
experiments. Numerical calculations with finer fc-meshes will be necessary to investigate 
the detailed momentum dependence of the order parameter on each small Fermi circle. 

We turn our attention to the anomalous Green's functions in the orbital representation, 
F£^£.^{k), which are obtained from Fa{k) by eq. (25). We display only three diagonal compo- 
nents in the orbital representation, Fxz xz{k), FyzYzik), and Fx2_y'^ x^-Y^ik)^ = ttT, 
since the other remaining diagonal elements are much smaller. Note that not all of the off- 
diagonal components in the orbital representation are so small, for example, FxzYz{k) and 
Fyz xz{k) (whose plots are not shown in the present article) have magnitudes comparable to 
Fxz xz{k) and FyzYzik). This is due to the fact that these orbitals are well hybridized with 
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(a) band 2 (d) XZ orbital 




(b) band 3 (e) YZ orbital 




Fig. 4. (Color online) Momentum dependence of the anomalous Green's functions at x = 0.10 and 
T =0.008 eV in the extended Brillouin zone, (a) F2{k,mT), (b) F3(fc,i7rr), (c) F4(fe,i7rr), (d) 
Fxzxz{k,inT), (e) FYZYzik,'nTT), and (f) F^^-y^ ^2_y2 (fc, IttT). Note that (a)-(c) showFa(fc) 
(a — 2,3,4) in the diagonalized-band representation, while (d)-(f) show three diagonal elements 
Fu{k) {£ — XZ, YZ, — Y'^) in the orbital representation. The indices a, a\ (3, 7, and 7' in (a)- 
(c) denote the Fermi pockets shown in Fig. 2(a). The green solid lines on the base plane indicate 
the zero contours of the anomalous Green's functions. We find in (a)-(c) that the green lines do 
not cross the Fermi surfaces, i.e., there are no nodes of the superconducting order parameter on 
the Fermi surfaces. 



each other. Fxz x'^-V^i^)-, z -V^ {^) 1 ^x'^-V^ xzi^)^ ^^^d -Fx2-y2yz(^) (whose plots are 
not shown in the present article) are negligibly small. An important feature of the anoma- 
lous Green's function in the orbital representation is that, roughly speaking, the maximum 
and minimum of Fu{k) in momentum space have nearly the same absolute magnitude but 
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Fig. 5. (Color online) Phase diagram. The horizontal and vertical axes denote electron doping level 
X and temperature (in units of eV). Negative values of x indicate hole doping. Circles (crosses) 
denote that the eigenvalue of the Eliashberg equation is larger (smaller) than unity, i.e, Amax > 
1 (Amax < !)• In other words, the circles and crosses represent superconducting and normal 
states, respectively. The solid curve represents the transition temperature determined by linear 
interpolation. The Coulomb integrals are U = 1.2 eV, U' — 0.9 eV, and J — J' — 0.15 eV. 



opposite signs to each other, and consequently Y2k-^uik), which indicates the strength of the 
on-site pair correlation at orbital I, is suppressed to a small value for each i. We stress that 
this is not an accidental feature but is essential for the unconventional s±-wave pairing due 
to electron correlations, as we shall discuss later in §5. 

3. 3 Doping dependence of transition temperature 

Experiments suggest that carrier doping is an important method for controlling iron- 
pnictide superconductivity. We investigate the effects of carrier doping on the transition tem- 
perature and pairing symmetry. In the present study, we assume simply that carrier doping, 
which is realized by the chemical substitution of atoms in experiments, only shifts the chemical 
potential and retains rigidly the band structure. The dependence of the transition temperature 
on doping level x is displayed in Fig. 5, where the transition temperature is determined by the 
linear interpolation of eigenvalues. We find that the superconducting state also extends in the 
hole-doped region x < 0. Thus, our theory is consistent with the existence of superconductiv- 
ity in hole-doped cases, which was suggested by some experiments. Although in practice the 
structural transition and SDW ordering inhibit the superconductivity near the undoped case 
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Fig. 6. (Color online) Doping dependence of the eigenvalues of the Eliashberg equation for U = 1.2 
eV, U' = 0.9 eV and J = J' = 0.15 eV at T 0.008 eV. 



X = 0, the structural transitions and the SDW ordering are unfortunately beyond the scope 
of our theory. 

Concerning the pairing symmetry, the spin-singlet s±-wave symmetry does not change 
over the whole doping region that we have calculated. The eigenvalues at T = 0.008 K for 
various symmetries are plotted in Fig. 6. The s±-wave state is likely to occur robustly against 
other pairing symmetries in iron pnictides. Varying the doping level, the volumes of the Fermi 
pockets are changed, but the zero contours of the anomalous functions always avoid crossing 
the Fermi surfaces. Thus, the fully gapped state is realized over the whole superconducting 
region. Here we note that the transition temperature is suddenly suppressed at a high doping 
level on both sides of the electron- and hole-doping regions. This is because the hole pocket 
a' around (vr, vr) disappears near this electron-doping level and the electron pockets 7 and 7' 
around (vr, 0) and (0, vr) disappear near this hole-doping level. 

3.4 Effect of vertex corrections 

It will be interesting to study the effects of the vertex corrections on the eigenvalues of 
the Eliashberg equation. The eigenvalues are compared for the following three cases: (i) the 
third-order calculation including the vertex corrections, using eq. (14), (ii) the third-order cal- 
culation without the vertex corrections, i.e., using A^^j^^C^) = ^CiL^^-' "''^CiL^^-' ~'~ ^CiC2^^^^^ 
instead of eq. (14), and (iii) the second-order calculation, i.e., using A(^-^^(^^(k) = A^^^^_^{k) + 
(2) 

CiC2^^^ instead of eq. (14). The calculated results are displayed in Fig. 7. We see that the 
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Fig. 7. (Color online) Eigenvalues of the Eliashberg equation for three cases: (i) third-order calcula- 
tion (including the vertex corrections), (ii) third-order calculation without the vertex corrections, 
and (iii) second-order calculation. The Coulomb integrals are U — 1.2 eV, U' = 0.9 eV and 
J = J' = 0.15 eV. 



second-order theory cannot quantitatively explain the high Tc within the realistic strength 
of the Coulomb interaction, because the eigenvalues are too small. The vertex corrections 
reduce the eigenvalues, and suppress the transition temperature. Thus, we consider that tran- 
sition temperatures are overestimated quantitatively, when using only the terms included in 
the RPA. This is qualitatively the same as in the d-wave superconductivity of single-band 
systems such as cuprates and organics.'^'^"^^^ The momentum dependence of the anomalous 
Green's functions Fa{k) calculated without the vertex corrections (whose plots are not shown 
in the present article) is almost the same as that calculated with the vertex corrections. 
We consider that the vertex corrections do not affect the momentum dependence of Fa(k) 
qualitatively, although they lower the transition temperature. 

3.5 Dependence on Coulomb integrals 

Eigenvalues of the Eliashberg equation for various sets of Coulomb integrals are displayed 
in Table I. We find that the s-wave (s±-wave) state is the most favorable pairing state robustly 
against the other pairing states. The order of pairing strength is s± > d^2_^y2 > p > dxy for 
all the cases of Coulomb integrals that we have calculated. In Table I, comparing the upper 
three cases of C/ = 1.2, we speculate that eigenvalues become small for small U' . We infer from 
this that processes including the interorbital Coulomb scattering U' play a role of enhancing 
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Table I. Eigenvalues of the Eliashberg equation for various Coulomb integrals, x = 0.10 and T 
0.004 cV. 



Coulomb integrals (cV) 



s-wave p-wave dx2_y2-wave dxy-wave 



u = 


1.2, U' 


= 0.9, J = 


J' 


= 0.15 


1.17 


0.45 


0.89 


0.14 


u = 


1.2, U' 


= 0.6, J = 


J' 


= 0.3 


0.56 


0.068 


0.39 


0.052 


u = 


1.2, U' 


= 0.9, J = 


J' 


= 


1.34 


0.58 


1.05 


0.18 


u = 


1.5, U' 


= 0.5, J = 


J' 


= 0.5 


1.09 


0.19 


0.71 


0.12 



the eigenvalues. By comparison with the last line in Table I, we find that the eigenvalues are 
reduced by using a small U', but are recovered by taking slightly larger U. 

4. Analysis of 2-band Hubbard model 

4-1 2-band tight-binding model by Raghu et al. 

We now use another simpler model. Raghu et al. proposed a two-band model, taking 
account of only the xz and yz orbitals:^^^ 

Exz{k) = 2ti cos kx + 2t2 cos ky -\- 4ts cos kx cos ky, (28) 

Eyz{k) = 2t2 cos kx + 2ti cos ky -\- 4t3 cos kx cos ky, (29) 

-E'hyb(fc) = — 4:t4 sin kx sin ky, (30) 

where the momentum variables kx and ky represent positions in the extended Brillouin zone, 
and the orbitals xz and yz should not be confused with the XZ and YZ orbitals in the original 
representation. £^hyb(^) is the hybridization term between the xz and yz bands. We take 
ti = 0.3, t2 = —0.39, t3 = 0.255, and t4 = —0.255, in units of eV, which give a quantitatively 
realistic energy scale, specifically, a total bandwidth of about 4 eV. The chemical potential is 
fi = 0.435 eV. The diagonalized band energies are 



Ei{k) = E+{k) - jE_{kf + E^y^{k)\ (31) 



E2{k) = E+{k) + ^JE^{k)^ + E^y^ik)^, (32) 

with 

Em = ^"'''t ''''' 
The electronic structure of this model is given in Fig. 8. The Fermi surface is qualitatively 
obtained from that by electronic band structure calculations and ARPES experiments, i.e., 
electron pockets around (0,0) and (7r,7r) and hole pockets around (7r,0) and (0,7r). How- 
ever, the energy band dispersions and the density of states are markedly different from those 
obtained from band calculations. 
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Fig. 8. (Color online) Electronic structure of the 2-band model by Raghu et al. in ref. 51. (a) Fermi 
surface, (b) energy band dispersion, and (c) density of states. 



4-2 Eigenvalues of Eliashberg equation 

For this simple 2-band model, we perform the same analysis as that for the 5-band model 
in §3. The calculated eigenvalues of the Eliashberg equation are presented in Fig. 9. The most 
favorable pairing symmetry is the odd-parity spin-triplet p-wave symmetry. This is due to the 
effect of the vertex corrections. As discussed in the case of Sr2Ru04,'^^^ the influence of the 
vertex corrections will be strong for such large particle-hole asymmetry as that in the present 
electronic structure. However the obtained eigenvalues are too small to explain the high Tc, 
and the triplet pairing state is not expected to be realized. Thus, we infer that the 2-band 
model is insufficient to describe the high-Tc superconductivity in iron pnictides. 

5. Discussion 

We have seen in §3.2 that the eigenvalues of the Eliashberg equation for the 5-band 
model are sufficiently large to explain the high transition temperatures. The highest transition 
temperature that can be expected according to Fig. 5 is about 160 K at around x = —0.1, 
and higher than the observed highest value of 56 K. It is, of course, possible to obtain more 
realistic values of Tc by taking slightly smaller values of the Coulomb integrals. In the present 
calculations, we have not included the normal self-energy corrections in the Green's functions. 
If we include these corrections, then Tc is expected to be somewhat decreased owing to the 
effect of quasi-particle damping and mass renormalization. For a more quantitative analysis, it 
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Fig. 9. (Color online) Eigenvalues of Eliashberg equation for the 2-band model as a function of 
temperature. The Coulomb integrals are U — 1.2 eV, U' — 0.9 eV, and J ~ J' ~ 0.15 eV. 



is necessary to include the normal self-energy corrections. In addition, it may be necessary to 
use more realistic three-dimensional electronic structures, since the observed resistivity does 
not show such strong two-dimensionality-'^^'^^-* as in other unconventional superconductors, 
such as cuprates, Sr2Ru04, and organics. On the other hand, the eigenvalues calculated for 
the 2-band model appears to be too small to explain the high T^. The 2-band model possibly 
lacks some essential factors, e.g., the xy states (equivalently the — Y^ states). In fact, as we 
have seen in the calculations using the 5-band model, the superconducting order parameter 
takes the maximum value on the X"^ — Y"^ state, and the inclusion of at least this orbital 
appears to be necessary to describe iron-pnictide superconductivity. 

We discuss the momentum dependences of the anomalous Green's functions. As we have 
seen in §3.2, the anomalous function changes its sign between the hole and electron pockets. 
Note that this is a natural result of electron correlations. In strongly correlated electron 
systems, the motion of electrons is generally affected strongly by the Coulomb repulsion. 
Electrons move in metals, avoiding each other, in order to reduce the expectation value of 
the Coulomb interaction. The most significant part of the local Coulomb interaction in actual 
solids is the on-site intraorbital interaction, i.e., the first term of the right-hand side of eq. (4): 

^'U = 2^ 4a4<x'Ci£a'Cifo- (34) 

In superconductors, the contribution of the superconducting condensate to the expectation 
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value of H'jj is evaluated within a simple mean-field approximation: 

SEu = {Hi) (35) 

= ^ i4ea4ea'(^ieT'Ciea) (36) 

^ (37) 

fefe' t o- 

where 

Fla<y{k) = {CkiaC-kla) (38) 

^L(fe) = (39) 

and Ckia is the Fourier transformation of Cj^^^. In odd-parity spin-triplet pairing states, SEu 
vanishes completely. We now consider even-parity spin-singlet pairing states. The spin part 
is decoupled by F£rra'{k) = F({k)[iay]a(T' and F^^^,{k) = {k)[—iay]aa' , and the expectation 
value of the energy is reduced to 

'^^^=^EiE^^(^)r- (40) 

e k 

It will be easy to see that, in order to reduce SEu, \ Ylk^eik)] should be small for each i. 
Here note that F({k) is calculated by the anomalous Green's functions:^^) 

Fe{k) = Tj2Fee{k). (41) 

The momentum dependence and nodal positions of Fu{k) (except absolute magnitudes) are 
generally almost the same at all frequencies. Therefore, we may expect that the momentum 
dependence of -F^(fe) is approximately the same as that of Fu{k). As we have seen in §3.2, 
the maximum and minimum of Fu{k) in momentum space have almost the same absolute 
magnitude but opposite signs to each other, and consequently Ylk^eeik) is suppressed to 
a small value for each £. Thus, we find that the calculated momentum dependences of the 
anomalous functions F£i{k) are indeed such that the on-site intraorbital correlation energy 
SEu is reduced. Using eqs. (25) and (27), we obtain the band representation 



(42) 



Y,T.\^ea{k)fFa{k). 

k a 

In many other unconventional superconductors including d3.2_j^2-wave superconductors such 
as some heavy-fermion superconductors, this value will be suppressed by the sign change of 
Fa{k) on each band a. On the other hand, in iron-pnictide superconductors, it is suppressed 
by the sign change of Fa{k) not on each band but among different bands. Thus, we naturally 
expect that the superconducting order parameter should always take nearly the same absolute 
value but opposite signs between the hole pockets around the F point and the electron pockets 
around the M point, in order to reduce the correlation energy SEu- The superconducting gap 
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function will also always take nearly the same absolute values between the hole pockets and 
the electron pockets, although the gap magnitude itself will depend on the temperature and 
doping level. This is expected to be confirmed by experiments, such as comprehensive ARPES 
studies on the band dependence of gap magnitude at various temperatures and doping levels 
in various iron-pnictide superconductors. If experiments indeed confirm this, they will support 
the realization of the s±-wave pairing symmetry in iron-pnictide superconductors. 

It will be interesting to compare our results with those of other theoretical studies. The sj-- 
wave or extended s-wave pairing states have been obtained by many other theoretical studies. 
Among them, the RPA calculations were performed by Kuroki et al. for the same 5-band 
Hubbard model^^^ (The 5-band model used in our present study was constructed by them), 
and by Yanagi et al. for a more complex d-p model including the As4j9 orbitals explicitly. 
We speculate that the s±-wave pairing state will be more robustly favorable against the next 
most favorable pairing state, i.e., the d^2_j^2-wave state, in the present perturbative study 
than in the RPA studies. For a higher doping level x = 0.3 (n = 6.3), where the hole pocket 
a' around (vr, vr) does not exist, the d,j.2_y2-wave state might become the most favorable as 
mentioned in ref. 23, but the eigenvalues will not be large enough to explain the high Tc- Our 
perturbative results suggest that the superconducting order parameter (and superconducting 
gap) does not have any nodes on the Fermi surface in the whole doping region, while nodes 
appear on the electron pockets around the M point in the RPA results of ref. 23, despite the 
fact that we have used the same model as Kuroki et al. used. The origin of this discrepancy 
is luiclear. Other more recent RPA results by Ikeda suggest that there are not nodes on the 
Fermi surface. Yanagi et al. suggested that the s-|--wave state is the most favorable in a 
d-p model with reasonable parameters, but the d^2_y2-wave pairing state (the dxy-wave state 
in the original folded Brillouin zone) becomes the most favorable, overcoming the s-t-wave 
state for small Hund's coupling or J = 0.^^^ On the other hand, in our present perturbative 
results using the Hubbard model, we could not find any sets of Coulomb integrals for which 
the dx2-y2-wave pairing state becomes the most favorable. We speculate that this difference 
originates from the difference in electronic structures used for the calculations (i.e., the 5- 
band Hubbard model or the more complex d-p model), rather than that in approximation 
methods (i.e., third-order perturbation theory or RPA). In fact, according to recent RPA 
results for the 5-band Hubbard model,^^) the dx2_y2-w&ve pairing state does not overcome the 
s-|--wave state even for J = 0, as obtained from our present third-order perturbation theory. 
Concerning the superconducting mechanism, we consider that the magnetic fluctuations due 
to the Fermi surface nesting between the hole and electron pockets are an important key 
factor, as considered in many other theoretical studies.^^'^^"^^) 

We can now understand why the transition temperature drops suddenly around the high- 
doping regions where the hole pocket around (tt, tt) or the electron pocket around the M 
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point disappears. This can be regarded as a natural result of the s±-wavc pairing due to 
electron correlations. For the higher doping levels, the Fermi surface nesting between the 
hole pockets and electron pockets cannot occur, because some of the nesting Fermi surfaces 
disappear. In addition, the sign change of the order parameter between them can no longer 
occur, and enhancement of the on-site correlation energy 6Ejj will inevitably destabilize the 
s-|--wave pairing state. We expect that the disappearance of superconductivity at the doping 
levels where the hole pockets or electron pockets disappear will be observed by systematic 
experiments in which both the Fermi surface and transition temperatures are measured for 
various doping levels. 

Finally, we mention that the 5-band model used in the present study is unsatisfactory 
regarding the following point. The volumes of the Fermi pockets appear to be quantitatively 
different from those of actual Fermi surfaces observed in ARPES experiments: the hole pockets 
have smaller volumes, while the electron pockets have larger volumes, than those in the ARPES 
results. This is the reason why the calculated transition temperature drops suddenly at a 
relatively low electron-doping level (about x = 0.18) and survives even up to a high hole- 
doping level (about x = —0.3), which is quantitatively inconsistent with experimental results. 
In actual systems, the superconductivity survives above the electron-doping level of x = 0.18, 
and vanishes at the high hole-doping level of x = —0.3. Thus, models possessing precise 
electronic structures and fermiology that are completely consistent with actual systems will 
be desirable for more quantitative studies. 

6. Conclusions 

In the present work, we have studied the high-Tc superconductivity in iron pnictides by 
using multi-band Hubbard models for iron 3d electrons and by solving the Eliashberg equation 
within the third-order perturbation theory. The unconventional s±-wave pairing symmetry is 
the most favorable and is highly likely to be realized, since it explains the high transition 
temperatures of iron pnictides. The superconducting gap is not expected to have any nodes 
on the Fermi surfaces in the whole doping region. This unconventional pairing symmetry 
is consistent with the results of many other theoretical and experimental studies. The mi- 
croscopic origin of this unconventional pairing is the local Coulomb interaction among the 
iron 3d electrons. Thus, we can regard this high-Tc superconductivity of iron pnictides as 
a natural result of the electron correlation, as in the cases of other established correlation- 
induced (non-phonon-mediated) superconductors, such as cuprate, ruthenate, organic, and 
some heavy-fermion superconductors. An important difference from other unconventional su- 
perconductors is that the order parameter changes its sign not on the Fermi surface but 
between the Fermi pockets. We hope that further comprehensive measurements such as those 
mentioned in §5 will be carried out in order to confirm the realization of this intriguing pairing 
symmetry in iron-pnictide superconductors. 
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